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TECHNICAL MEMORANDUM 78141 


A TECHNIQUE FOR SIMULATING TURBULENCE FOR AEROSPACE 
VEHICLE FLIGHT SIMULATION STUDIES 

INTRODUCTION 

The purpose of this document is to describe a technique for simulating 
atmospheric turbulence which involves a large variability of turbulence proper- 
ties durit^g the simulation procedure. The technique, which has direct applica- 
tion to Si.attle Orbiter entry and aircraft landing simulation studies, involves 
extension of recursive filter techniques which have been widely used in the 
aerospace and aeronautical engineering communities. To apply these filters, 
a white noise process is passed through a linear filter to generate output 
processes with statistics that resemble those of atmospheric turbulence (1-31. 

In most applications the simulated turbulence velocity components are Gaussian 
processes with one -dimensional spectra with Dryden or von Karman mathe- 
matical forms (4,5) . The defining parameters of these spectra are the 
turbulence velocity standard deviations and integral scales. These parameters 
also are the defining parameters of the filters used to process the white noise 
to generate the turbulence. In fact, the filters are manufactured from the 
spectral shapes of the turbulence velocity spectra (3| . In many applications the 
standard deviations and integral scales of turbulence can vary significantly during 
a flight simulation. For example, from Shuttle Orbiter entry at 120 000 m down 
to 10 000 m the standard deviations and integral scales of the longitudinal, 
lateral, and vertical components of turbulence vary by one to two orders of 
magnitude [6J. Furthermore, from 10 000 m down to 18.3 m these parameters 
vary by one-half to two orders of magnitude, with the largest variation taking 
place below 300 m. To apply the linear filters in this situation, the turbulence 
velocity integral scales and standard deviations are varied during the flight 
simulation in accordance with a priori specified functions of altitude. Thus, 
during a flight simulation, the flight trajectory and the turbulence simulation 
procedure are tied together in the sense that the flight trajectory defines altitude 
which serves to define the turbulence integral scales and standard deviations 
which determiiK the simulated turbulence, which, in tura, feeds back to the 
vehicle and affects the trajectory. Thus, to accommodate the relationship 
between the flight trajectory and the turbulence simulation model, the turbulence 
must be simulated during the flight simulation. This can place large demands 
on the computer facilities available for the flight simulation and could preclude 
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a high fidelity flight simulation. However, this problem can be circumvented 
by using simulated turbulence generated prior to the actual flight simulation. 

This can be accomplished via two procedures. One procedure involves the use 
of a ’’nominal” trajectory to define the variability of the turbulence standard 
deviations and integral scales during the flight simulation. The second procedure, 
the one with which this report is concerned, involves the use of preflight simu- 
lated turbulence time histories which are in nondimensional format wherein the 
turbulence velocity components are scaled with the appropriate standard devia- 
tions and the time coordinate is scaled with the appropriate integral scale 
divided by the true air speed of the vehicle. The central idea here is that the 
nondimensional txmbulence is applicable anywhere along the trajectory. Thus, 
during the flight simulation, the nondimensional turbulence is placed into the 
flight simulation computer (via data tape, for example), and the dimensional 
turbulence as a fimction of dimensional time is generated by applying the scaling 
factors (integral scale and standard deviation) and true air speed as a function 
of time along the trajectory. Application of the scaling factors ties the simu- 
lated dimensional turbulence to the flight trajectory with a potentially large 
saving of flight simulation computer usage. The purpose of this report is to 
describe this nondimensional turbulence simulation concept and to discuss its 
application to aerospace and aeronautical flight simulation. 


TURBULENCE MODEL 


The turbulence simulation technique described in this report is based on 
a Gaussian turbulence model in which the spectra for the longitudinal, lateral, 
and vertical components of turbulence are given by 


L a" 

Longitudinal: (a,) = f 1 ^(l7m/V) ’ 

L 2a* 1 + 3(L w/V)* 

Lateral: M ‘ TT; (1 J7v)’|' <2) 


Vertical: (f>^ (co) 


L 2a 2 1 + 3(L o)/V)* 

w ^ w ' 

V >r [1+ (I^w/V)*]* 


( 3 ) 
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where L and a denote the integral scale and standard deviation of turbulence, 
CO is tlie radian frequency, and V is the true air speed of the vehicle. The 
reader is referred to Section 5. 3. 14 of Reference 5 for further details. The 
quantities L , o , and V are functions of altitude and, hence, time for a non- 
horizontal flight path; however, they will be treated as being locally constant 
along the trajectory to generate the turbulence filters and then will be permitted 
to vary along the trajectory to accommodate the nonstationary character of 
turbulence resulting from a vehicle traversing a nonhomogeneous turbulence 
field. 


TRANSFORMATIONS 


By transforming ctJ to a nondime ns ional frequency 


Q 


Lct> 

V 


( 4 ) 


and transforming the dimensional spectral density function (co) to a non- 
dimensional spectral density 


$(n) 



<P 





( 5 ) 


equations ( 1) , ( 2) , and ( 3) can be written as 


^u^^u^ " TT 1 + n * 
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( 6 ) 
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V”v^ = 7 li + n *1* 
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where the subscripts on fl refer to the fact that L it L # L . The afore- 

u V w 

mentioned spectra imply the existence of Fourier decompositions of the form 


“a ^ 

u(g = / u(ge "dn^ 

-iCO 

» A ‘«v'v 

V(g = / V(S!^e 

® A f 

W(S ) = / W(n ) e "^dfl 
' w' ' w w 

-00 


( 9 ) 

( 10 ) 

( 11 ) 


where (U) denotes Fourier decomposition of (U) over the domain 
The quantity f is a nondimensional coordinate, namely 


t 


tv 

L 


(12) 


where t is dimensional time measimed with respect to an arbitrarily selected 
point on the vehicle entry trajectory. The quantities U , V , and W at the 
nondimensional longitudinal, lateral, and vertical components of turbulence are 
defined as 


«(t„) 


u(t L /V) 
' u u ' 


or 
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( 13 ) 


V(t ) 
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'(t, 


( 14 ) 
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( 15 ) 


W(£ ) = 
w' 


*(«w 
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where u , v , and w are the dimensional longitudinal, lateral, and vertical 
components of turbulence which are functions of t . Note that after applying the 
transformation, equation (12), we have expressed nondimensional components 
of turbulence U , V , W as functions of nondimensional time £ measured with 
respect to an arbitrarily selected point on the Orbiter t. ijectory. In view of the 
previous developments, it appears reasonable to simulate turbulence in terms 
of the nondimensional random functions U(£^) , V(f^) , and W(£^) because 

the variability in the turbulence simulation resulting from time dependent ct, L , 
and V has been absorbed into nondimensional random functions U(S^), V(5^), 

and W(£ ) so that a, L, and V do not appear explicitly in a nondimensional 

context. 


DIGITAL SIMULATION FILTERS 


To mechanize the concepts enunciated in the previous sections we shall 
use the concept of digital recursive filters. We shall use the digital filters 
associated with Dryden spectra, given in Reference 7. Thus, the components 
of the digitized turbulence are simulated by passing white noise through the 
following filters; 


U 

n 


C U , 
u n-1 


+ D 

u 


N 

u,n-l 


(16) 


V 

n 


(17) 


W = C. 




(18) 


where subscript n is a running index (i.e., n= 1, 2, . . . ), and N^, N^, and 
N^ are independent Gaussian white noise sources. The remaining quantities 
are defined as 
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Do = k ^ ^ 
2,w w a * 
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-b -a A£ 

— - A? e ^ ^ . (28) 
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The various parameters in equations (19) through (28) are defined as follows: 


k = - 

U 7T 


k = k 
V w 


and Af » Af , and A^ denote the nondimcnsional time intervals between 
u V w 

the digitized values of U, V, and W, respectively. 

The procedure for simulating turbulence consists of (1) selecting values 
of U, * V , V , W , and W (by random selection in the computer, for example) , 

(2) selecting values of N ,, N ,, N N ,, and N (3) substituting 
' ' ^ u,l V, 1 v,2 v,l v,2 ' ' 

results of steps (1) and (2) into equations (16), (17), and (18) to calculate 

U„» V„, W , (4) incrementii^ the subscript n by 1 and repeating steps (1) 

2 u u 

through (3) a sufficient number of times to generate the desired nondimensional 
time histories of turbulence, and (5) transforming to dimensional time histories 
with equations (12) through (15). 

The digital simulation of U, V, and W requires specification of the non- 

dimensional intervals A£ , A£ , and At and, of course, three independent 

Ul V IV 

Gaussian white noise processes. The dimensional interval At can be written as 


At = ^ AJ . 
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Furthermore, the interval At can be written in terms of the local Nyquist 
freqtinncy a? of the data set, 



(J4) 

# 


Combining equations (33) and (34) yields 


O). 


N 



where 


4< = 


nV 

L • 


(35) 


(36) 


During a nondimensional turbulence simulation is held constant so that 
according to equation ( 35) the Nyquist frequency will vary in direct proportion 
to il > . Thus, to encompass all turbulence frequencies of importance during a 
flight simulation we specify Af by selecting the largest value of of concern 
and the smallest value of ip, (w^^) • and j/^. Thus, 


!/■ 


= 


(oj ) 

^ n'o 


(37) 


APPLICATION TO ORBITER ENTRY 



As stated in the previous section, the simulation technique requires 
specification of At . Let us consider the entry of the Shuttle Orbitcr as an 
example* Figure 1 contains a plot of the Shuttle Orbiter true air speed V during 
the entry flight phase as a function of altitude above mean sea level (MSL) for 
the first Operational Flight Test (OFT) . Figures 2 and 3 contain plots of the 
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Figure 3. Shuttle program natural environment design standard deviations of 
turbulence (a) as a function of altitude [6) . 



integral scale of turbulence (L) and the standard deviation (a) for each com- 
ponent of turbulence as specified by the Shuttle Level II requirements document 
[6], As stated earlier, a serves as velocity scale and L/V serves as a time 
scale to nondimensionalize the turbulence velocity time histories. Combination 
of the data in Figures 1 and 2 yields the calculated values of ip [equation (36) ] 
in Table 1 for each component of turbulence. As discussed in the previous 
section, we require the lowest value of ip and the largest value of w to 
encompass all turbulence frequencies of concern. 


TABLE 1. VALUES OF p DERIVED FOR OFT No. 1 MISSION 


Altitude (m) 

' ^ 

w 


4> 

^u 

120 000 

1. 884956 

0. 113097 

0.113097 

90 000 

2.748894 

0.169163 

0.169163 

70 000 

3.257948 

0. 204569 

0.204569 

50 000 

3.337942 

0.205412 

0.205412 

30 000 

2.094395 

0. 142800 

0.142800 

20 000 

1.636246 

0.261799 

0.261799 

18 000 

1.570796 

0.329119 

0.329119 

16 000 

1.374447 

0.407243 

0.407243 

14 000 

1.308997 

. 0.496041 

0.496041 

12 000 

1.393448 

0. 863938 

0.863938 

10 000 

1.481883 

1.481883 

1.481883 

7 000 

1.274420 

1. 274420 

1. 274420 

5 000 

1.185507 

1. 185507 

1.185507 

2 000 

1.037318 

1.037318 

1.037318 

1 000 

0.978043 

1.978043 

0.978043 

500 

1.2566 

1.2566 

1.2566 

300 

1.6232 

1.6232 

1.6232 

50 

12.5664 

8. 7965 

6. 1087 

18.3 

37.6991 

20. 9440 

11.7810 
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According to Table 1, the values of listed in Table 2 can be used for 

the Shuttle Orbiter mission. Table 1 shows that if) is very nearly equal to a 
constant for the vertical component of turbulence for the entry mission except 
in Ihe vicinity of the ground where 4 ' increases rapidly with decreasing altitude. 
However, for the longitudinal and lateral components of turbulence, 4> varies 
by one order of magnitude from 120 000 m down to 10 000 m and remains 
constant imtil, in the vicinity of the ground, ip increases rapidly with decreasing 
altitude. The actual Nyquist frequency for each component of turbulence in the 
data set will vary with altitude according to the formula 


a; = (c ) -^ 
N ^ N^o 


(38) 


so that any variability of ip with altitude will be reflected in Thus, for 

N 

the vertical component of turbulence, is very nearly constant with altitude 

except near the ground, where 4 and, hence, increase by nearly one order 

of magnitude. For the longitudinal and lateral components, the respective values 
of increase by one order of magnitude in the altitude band from 120 000 to 

10 000 m, as well as near the ground. 

TABLE 2. VALUES OF 4 > FOR TURBULENCE 

o 

SIMULATION DERIVED FOR THE ORBITER 
OFT No. 1 MISSION 


Component I 

*^0 

Longitudinal 

0.113 

Lateral 

0.113 

Vertical 

0.978 


A value of > 1 implies that the digital rate at which turbulence 

N ^ o 

is provided in the simulation is greater than the rate which would actually be 

necessary. Thus, if data tapes are used to store the digital simulations of 

turbulence, values of 41/41 > 1 and, hence, a)^/(co ) > 1 near the ground 

O N N n 
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do not appear to pose data storage problems for the vertical component of 

turbulence because the length of the trajectory near the ground, over which the 

data rate would be too large, is relatively short. However, in the case of the 

longitudinal and lateral components of turbulence, varies by an order of 

N 

magnitudt) . ver a significantly large portion of the entry mission, i.e. , 120 000 
to 10 000 m (Table 1). Thus, 0,113 corresponds to a digital simulation 

rate ahiO . is too large by a factor of 10 below 10 000 m down to 50 m and by a 

factor of 100 below 50 m. To remedy this situation two values of ^ could be 

used lor the longitudinal and lateral turbulence simulation (Table 3). The 

increase in $/il> , and hence w /(w ) , by an order of magnitude below 
O N N Q 

500 m dor s not appear to pose data storage problems for the same reason that 
the vertir al component does not ^pear to pose data storage problems relative 
to simi:i.'.' ed turbulence near the groimd, as previously stated. 


TABLE 3. VALUES OF ip^ FOR LONGITUDINAL AND LATERAL 

TURBULENCE SIMULATION AS A FUNCTION OF ALTITUDE 
DERIVED FOR THE ORBITER OFT No. 1 MISSION 


Altitude > 10 km 

Component 

^0 

Longitudinal 

0.113 

Lateral 

0.113 

Vertical 

0.978 

Altitude < 10 km 

Component 

^0 

Longitudinal 

0.978 

Lateral 

0.978 

Vertical 

0,978 


1-v 
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TURBULENCE DATA TAPE CONCEPT 


To provide nondiinensional turbulence for flight simulation of the Shuttle 
Orbiter, a data t^e with nondimensional turbulence for at least two sets of 

values would be required. First, an upper bound Nyquist frequency (cOj^) would 

N o 

be required. This upper bound frequency should be a factor of 10 greater than 
any turbulence frequency of concern, in an engineering context, to insure 

that all Fourier components of the turbulence with frequency w < w have 

100 percent (or very nearly so) spectral energy. The spectral densities <p (w) 

s 

of the simulated turbulence for w > will tend to zero as cu — (w.,) and 

t No 

are related to the spectral densities of the turbulence model [equations (l) , (2), 
and (3)1 ?>(w) as follows: 




V^^E/ 


0(w) / itu 


K) 


Once the upper bond Nyquist frequency (co ) and ip are specified, the non- 

N o ® 

dimensional time intervals [equation (37) ] are known and equations (16), 

( 17) , and ( 18) can be applied. 

To provide nondimensional turbulence data tapes to encompass the total 
Shuttle Orbiter entry mission, the tapes could consist of two parts in sequence, 
or more parts if required for more detailed definition of the turbulence with 
respect to altitude. The first set of nondimensional simulated turbulence 
velocities is applicable to altitudes above 10 000 m and is generated with the 
values of ip^ given in Table 2. The second set is applicable to altitudes below 

10 000 m and is generated with ip^ - 0.978 for each turbulence component. 

Thus, during a flight simulation, nondimensional turbulence from the first 
turbulence set on the data tapes would be used down to 10 000 m whereupon the 
data tape could provide for a number of different simulations by entering sets 
1 and 2 at different locations with the actual number of different simulations 
available depending on the relative magnitudes of ip^ and (Wj^) . 








In each data tape set the nondimensional time between the simulated non- 
dimensional turbulence velocities for each component of turbulence is a constant. 
In view of equation (33), the dimensional time At between the dimensional 
simulated winds will vary with altitude to the degree that 1,/V varies with 
altitude. Furthermore, At for the various components of turbulence will not 
necessarily be equal. Above 10 000 m ( At)^ = (At)^ # ( At)^. Below 10 000 m 

and above 304.8 m (At)^= (At)^= 304.8 m (At)^ # (At)^ # 

(At) . However, this variation of At from component to component and from 

altitude to altitude can be easily accommodated during the flight simulation 
process. 
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